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ABSTRACT: 

This paper presents an optimization of corridor observation method ( COM) which is an 
applicable optimization algorithm based on the evolutionary algorithm to solve an environmental and 
economic Dispatch (EED) problem. This problem is seen like a bi-objective optimization problem 
where fuel cost and gas emission are objectives. In this method, the optimal Pareto front is found using 
the concept of corridor observation and the best compromised solution is obtained by fuzzy logic. The 
optimization of this method consists to find best parameters (number of corridor, number of initial 
population and number of generation) which improve solution and reduce a computational time. The 
simulated results using power system with different numbers of generation units showed that the new 
parameters ameliorate the solution keep tier stability and reduce considerably the CPU time ( time is 
minimum divide by 4) comparatively at parameterization with originals parameters. 
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I. INTRODUCTION 

The economic dispatching (ED) is one of the key problems in power system operation and planning. 
The basic objective of economic dispatch is to schedule the committed generating unit outputs so as to meet the 
load demand at minimum operating cost, while satisfying all equality and inequality constraints. This makes the 
ED problem a large - scale highly constrained non-linear optimization problem. In addition, the increasing 
public awareness of the environmental protection and the passage of the Clean Air Act Amendments of 1990 
have forced the utilities to modify their design or operational strategies to reduce pollution and atmospheric 
emissions of the thermal power plant. 

Several strategies to reduce the atmospheric emissions have been proposed and discussed. These 
include: installation of pollutant cleaning equipment, switching to low emission fuels, replacement of the aged 
fuel-burners with cleaner ones, and emission dispatching. The first three options require installation of new 
equipment and/or modification of the existing ones that involve considerable capital outlay and, hence, they can 
be considered as long-term options. The emission dispatching option is an attractive short-term alternative in 
which the emission in addition to the fuel cost objective is to be minimized. Thus, the ED problem can be 
handled as a multi-objective optimization problem with non-commensurable and contradictory objectives. In 
recent years, this option has received much attention [1-5] since it requires only small modification of the basic 
ED to include emissions. 

In the literature concerning environmental/economic dispatch (EED) problem, different technics have 
been applied to solve EED problem. In [1, 2] the problem was reduced to a single objective problem by treating 
the emission as a constraint. This formulation, however, has a severe difficulty in getting the trade-off relations 
between cost and emission. Alternatively, minimizing the emission has been handled as another objective in 
addition to the cost [5]. However, many mathematical assumptions have to be given to simplify the problem. 
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Furthermore, this approach does not give any information regarding the trade-offs involved. In other research 
direction, the multi -objective EED problem was converted to a single objective problem by linear combination 
of different objectives as a weighted sum [3], [6]. The important aspect of this weighted sum method is that a set 
of non-inferior (or Pareto-optimal) solutions can be obtained by varying the weights. Unfortunately, this 
requires multiple runs as many times as the number of desired Pareto-optimal solutions. Furthermore, this 
method cannot be used in problems having a non-convex Pareto optimal front. To overcome it, certain method 
optimizes the most preferred objective and considers the other objectives as constraints bounded by some 
allowable levels [5]. The most obvious weaknesses of this approach are that, they are time-consuming and tend 
to observe weakly non-dominated solutions [5]. 

The other direction is to consider both objectives simultaneously as competing objectives. The recent 
review to the Unit Commitment and Methods for Solving [7] showed that evolutionary algorithms are the most 
used in this case; certainly because they can efficiently eliminate most of the difficulties of classical methods 
[5]. The major problems of these algorithms, is to find the Pareto optimal front, to conserve the non- dominated 
solutions during the search and relatively long time to find the solution. 

In [8] we have proposed one method, based on the evolutionary method where the optimal Pareto front 
is obtained by the concept of corridor observation and the loses of non-dominated solution is reduce by the 
dynamism of archives during the different generations. The quality of solution and CPU time depend to the 
number of corridors, the initial population and the number of generations. In this paper, we propose some 
parameters which keep stable the solution and reduce considerably CPU time of COM. In the second part of this 
paper, we present materials and methods to solve the EED problem and in the third part, simulation and results 
are presented to enable us to find the new parameters and demonstrate their effectiveness by comparing it with 
the original settings. 

II. MATERIALS AND METHODS [8] 

In this part, we formulate the EED problem and present our approach to solve it. 

2.1. Problem formulation 

The EED problem is to minimize two competing objective functions, fuel cost and emission, while satisfying 
several equality and inequality constraints. Generally the problem is formulated as follows: 

2.1.1. Problem objectives 

> Minimization of fuel cost 

The generator cost curves are represented generally by quadratic functions. The total fuel cost ($/h) in terms of 
period T, can be expressed as: 

T Ng 

F (Pu)= Z Z '/' )> + ST. t (l- /„,_,)/„, (1) 
,=i ,=i where 

C /,(P,.,) = a i + b ,P ,,, + C ,P^, 
(2) 

andc „ (p. t ) is the generator fuel cost function; a ; , b and C; are the cost coefficients of i th generator; P it , is the 

electrical output of i th generator; Ng is the number of generators committed to the operating system ; Ij. t the statut 
of different generators; ST; the start-up cost . 

> Minimization of gas emission 

The atmospheric pollutants such as sulphur dioxides (SO2) and nitrogen oxides (NOx) caused by fossil-fuelled 
thermal units can be modelled separately. However, for comparison purposes, the total emission (ton/h) in one 
period T of these pollutants can be expressed as: 

T N S 

E{p,. l )=Y J Y j e fl (p il )I ij (3) 

r=l 1=1 

where 

e fl {P it ) = (a, + p.p tl + S iP l) (4) 
and a;, Pi, §i are the emission coefficients of the i lh generator. 

2.1.2. Objective constraints 

> Power balance constraint 
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Ng 



^ Spinning reserve constraint : 



I. < 0 



max, i i,t 



> Generation limit constraints : 

P,„„, , 7 ,,, ^ P,,t ~ p »«. >'<,, ' ' = 1 N S 

> Minimum up and down time constraint : 



(5) 

(6) 
(7) 



I i if t°" -< t; 

\ 0 if T° ff < T* 
0 or 1 otherwise 



(8) 



Where Ti Up represent the minimum up time of unit- i ; T^ 0 "" the minimum down time of unit i T° is the 

continuously off time of unit i and T; 0 " the continuously on time of unit- i . 
y Start-up cost 

f HST if T." < T" ff < T." + T""" 



ST 



I CST . if T° n y T. cau + T, ' 



2.2. The corridors observation method [8] 

The different steps of the COM are presented in the following figure 

/ Generate initial \ 
population J 



Evaluate function cost 
Evaluate emission cost 



Segment objectives space function 
in corridor 



Form a new population 



Search the best persons in 
corridors 



Archive the best persons 



Apply cross and mutation operator 




determine the pareto front 



determine the best compromize 
solution 



End 



Figure 1: Different steps of the algorithms [8] 

> Step 1 

In the first step, we start with to the status of different units generation, where we create randomly the initial 
population. Each individual is a combination of each power generation unit 



> Step 2 

In the second, using equations (1) and (3) we evaluate the objective functions of this population 
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> Step 3 

Using the minimum of the different objective functions of individuals who respect the constraints (5) to 
(8), we define the space solution and segment it to the corridors observations following the different axes which 
are specify by each function. 

> Step 4 

In each corridor, we search the best individuals who have the minimum objectives functions, and the 
non feasible solutions are classified using the number and the rate of violation constraints. Those solutions will 
be used to increase the number of feasible solutions. 



> Step 5 

We keep in the archives those best individuals 



> Step 6 

We verify the stopping criteria define as [10]: 

f = H d ) 

where 



(9) 





cl f 


1 






I 


CI \ 









(10) 

Explain the metric progression of the best individuals in each corridor. N F is the number of objectives 
functions ; Cl the number of corridor ; , F 4 ~ y the j" 1 objective function of the best individual in i th corridor ; 
F min and F max the minimum and maximum of the j function ; t is the present generation, t-1 the anterior 
generation. At times the maximum number of generation can be the alternative stopping criteria 



> step 7 

If the stopping criteria is not verified, we construct the new population using the selection, cross and 
mutation operators apply to the archive population and we return to step 2. 

> Step 8 

If the stopping criteria is verified we find the best compromise solution among the individuals of the 
Pareto front. Due to imprecise nature of the decision maker's judgment, each objective function of the i-th 

solution is represented by a membership /u . function defined as 

1 if F. < F.""° 
F™ -F 

p. = \ ■ — if F < F < F."" 

j F ma - p ' ' 

For each non-dominated solution, the normalized membership function is /u 1 calculated as : 



^ M N F 

-=i (12) 
where M is the number of non-dominated solutions. The best compromise solution is the one having the 

maximum of ju k . 



III. SIMULATION AND RESULTS 

In order to find the new parameters of COM and solve the EED, a 3 -units generation system is tested 
[11]. And extent at 6, 10 and 15 units generation. These parameters are applied in COM with the same software 
and computer used in [8]. The results are compared with the originals settings of COM. The data concerning 
the units generation are given in tables 1 and 2 [11]. 
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Table 1. The 3-units system data 








Unit 


p max 

(MW) 


p min 

(MW) 


a b 

($ /h) ($ /MWh) 


c 

($ /MW 2 h) 


R,"" 

($ /MW) 


down 

($ /MWh) 




1 


600 


150 


561 7.29 


0.00156 


100 


100 




2 


400 


100 


310 7.85 


0.00194 


80 


80 




3 


200 


50 


78 7.97 


0.00482 


50 


50 








Table 2. SO : 


and NOx coefficients emission gas data of 3-units 




Units 


a ; j (ton/h) 


(ton/h) 


P w 2 (ton/MWh) 


P Nox (ton/MWh) 


8 m 2 (ton/MW 2 h) 


(ton/MW 2 h) 


1 




0,5783298 


0,04373254 


0,00816466 


-9,4868099 e' 6 


l,6103e* 


1,4721848 e 7 


2 




0,3515338 


0,055821713 


0,00891174 


-9,7252878 e" 5 


5,4658 e 6 


3,0207577 e" 7 


3 




0,0884504 


0,027731524 


0,00903782 


-3,5373734 e" 4 


5,4658 e 6 


1,9338531 e 6 



In the implementation, we add the different coefficients of each gas per groups to have the coefficient of the 
whole gas 

3.1. Study of objectives according to the number of generation 

The study of convergence with the same parameters of [8] i.e., 50 corridors, 300 number of initial population is 
presented in figure 2 



cost convergence with generation 
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. 9270 
' 9265 
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'1 
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Figure 2: Convergence of fuel cost ($) and gas emissions (ton/h) objective functions (1000MW). 

These curves convergence show that since the 500 generations fuel cost and gas emission is uniform so, 
if we choose the number of generations at the value greater or equal than 500 generations the solution will be 
the same. So we can take 500 like the new number of generations. 

3.2. Study of objectives according to the number of corridors 

The representation of objectives, according to the number of corridor, with initial population number 
and number of generations set respectively at 300 and 1000 is shown in figure 3. 



www.ijceronline.com 



Open Access Journal 



Page 5 



Optimization of Corridor Observation Method to Solve Environmental and Economic ... 




50 60 
number ofcondois 



£.7525 




I 7 52 



number of castas 



Figure 3: Evolution of fuel cost ($) and gas emissions (ton/h) according to the number of corridors (850 MW). 
These curves shows that from 20 corridors fuel cost and gas emission are almost uniform. The Pareto front at 
this value of number of corridors is represented in figure 4. 
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Figure 4: Pareto front with the new parameters. 
Comparing with the Pareto front obtained with the original parameters [8] in figure 5 
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Figure 5: Pareto front with the original parameters. 
We can conclude that the increase of corridors number improves the Pareto front but from 20 corridors, the 
optimal solution is almost uniform. So we can consider 20 like new number of corridor. 



3.3. Study of objectives according to the number of initial population 

The representation of objectives, according to initial population, with number of corridor and number of 
generations set respectively at 50 and 1000 is shown in figure 
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Figure 6: Evolution of fuel cost ($) and gas emissions (ton/h) according to the number of initial individuals 

(850MW) 

These representations show that from an initial population of 200, fuel cost and gas emission are almost 
uniform. So we can set this parameter at 200. 

3.4. Comparative study between the original parameters and new 

To present the effectiveness of the new parameters in COM to unit commitment and EED, we applied it 
to a production plan of 3 -units during 5 hours and made the comparison study with the original setting. 
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Table 3. Unit commitment and EED during 5 hours with the two groups of parameters 



Hours Demands 


PI 


P2 


P3 


Fuel cost 


Emission 


Average 


(H) (MW) 


(MW) 


(MW) 


(MW) 


X10 3 ($) 


gaz 


CPU time 












(ton/h) 


(s) 


With original parameters 


1 550 


0 


0 


549.9 


5.0415 


5.6380 




2 600 


0 


0 


599.9 


5.4957 


6.1468 




3 650 


108.0382 


0 


541.8742 


5.9646 


6.6979 


23.29 


4 700 


114.5484 


0 


585.351 


6.4169 


7.2059 




5 750 


150.228 


0 


599.67 


6.8777 


7.7323 






With new parameters 


1 550 


0 


0 


549.9 


5.0417 


5.6383 




2 600 


0 


0 


599.9 


5.4957 


6.1468 




3 650 


100.5037 


0 


549.39 


5.9647 


6.6967 


5.1 


4 700 


106.2763 


0 


593.6432 


6.4179 


7.2054 




5 750 


152.6171 


0 


597.32 


6.8787 


7.7343 





The findings of this table are as follows: in terms of unit commitment, the results are identical; sensibly 
the same in terms of EED but the CPU average times is considerably reduced. This study have extended with 6, 
10, and 15 unit in table 4 findings was the same but the convergence speed according to the number of units is 
reduce with news parameters 



Table 4. Comparative study with 3, 6, 10 and 15 units generation 



Number of units 


3 


6 


10 


15 


With original parameters 


Fuel cost x 10 3 ($) 


9.2606 


9.1911 


9.176 


9.177 


Gas emission! ton/h) 


10.7347 


10.2771 


10.3679 


10.3683 


Average times (s) 


23.29 


31.92 


49.40 


104.61 


With new parameters 


Fuel cost xl0 3 ($) 


9.2583 


9.1922 


9.1844 


9.1806 


Gas emission (ton/h) 


10.7617 


10.1988 


10.2854 


10.3898 


Average time (s) 


5.14 


5.74 


7.03 


34.30 



IV. CONCLUSION 

In this paper, news parameters are proposed to optimize COM. COM have been proposed in [8] to 
solve environmental/economic power dispatch optimization problem and unit commitment. The study of 
objectives according to the number of generation, of corridors and initial population have allowed us to propose 
new parameters that are 20 corridors, 200 individuals for initial population and 500 generations. The 
parameterization of COM with these settings conserves the quality of solution in terms of unit commitment and 
EED. The principal advantages of this parameterization are the reduction of CPU time (the time is minimum 
divided by 4) and the convergence speed according to the number of generation units comparatively at the 
parameterization with originals settings (50 corridors, 200 individuals and 1000 generation). 
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